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Abstract 

We consider a reflected Ornstein-Uhlenbeck process driven by a fractional Brownian motion 
with Hurst parameter H £ (0,f). Our goal is to estimate an unknown drift parameter a £ 
(—oo, oo ) on the basis of continuous observation of the state process. We first derive two kinds 
of fractional Girsanov formula for H £ (0, 1) and show that they are equivalent, by establishing a 
transformation result between a fractional Brownian motion and a standard Brownian motion. 
Then, we investigate the standard maximum likelihood estimator of the drift parameter and 
prove its strong consistency and asymptotic normality. As an improved estimator, we derive the 
explicit formulas for the sequential maximum likelihood estimator and its mean squared error 
by assuming the process is observed until a certain information reaches a specified precision 
level. The estimator is shown to be unbiased, uniformly normally distributed, and efficient in 
the mean square error sense. 

Keywords: Reflected fractional Ornstein-Uhlenbeck processes, fractional Brownian motion, frac- 
tional calculus, parameter estimation, maximum likelihood estimator, sequential maximum likeli- 
hood estimator. 

AMS Subject Classifications: Primary 60G22; secondary 60H30, 62M09, 90B22 

1 Introduction 



We consider a drift parameter estimation problem for a one-dimensional reflected fractional Ornstein- 
Uhlenbeck (RFOU) process with infinitesimal drift —ax and infinitesimal variance a 1 , where 
a £ (—00,00) and a > 0. The RFOU process can serve as approximating models in diverse 
applications such as in physical, biological, and mathematical finance models (see, e.g., [25| [Tj HI [5] 
and also Section [2] below). The RFOU behaves as a standard FOU process in the interior of its 
domain (6,00). However, when it reaches its boundary b, then the sample path returns to the 
interior in a manner exercising with minimal "pushing" force. Our main interest in this model 

* E-mail: chihoon@stat .colostate.edu 

* E-mail: txjsong@hku.hk 



1 



stems from the fact that the RFOU process arises as the key approximating process for queueing 
systems with reneging or balking customers with long range dependent inter-arrival and/or service 
time processes (see [27\ W$\ [29] and the references therein). In such cases, the drift parameter a 
carries the physical meaning of customers' reneging (or, balking) rate from the system. More details 
are provided in Section [2] with regards to how the RFOU model can arise in the applications. In 
this paper, we expand the previously known results on the parameter estimation problems for a 
reflected Ornstein-Uhlenbeck (ROU) process to the case when the noise process is given by a frac- 
tional Brownian motion with Hurst parameter H € (0, 1). Such results require several nontrivial 
technical efforts, since the fractional Brownian motion is not a Markov process nor a semimartingale 
(unless H = 1/2) and the classical stochastic calculus is inapplicable in its analysis. 

We describe the RFOU model more precisely. Let A := (fi, J 7 , (J r t)t>o ) P) be a complete prob- 
ability space with the filtration (J-t)t>o satisfying the usual conditions. Define the RFOU process 
{Xt : t > 0} reflected at the boundary b £ M + on A as follows. Let {Xt : t > 0} be the strong 
solution (whose existence is guaranteed by an extension of the results of |17j ) to the stochastic 
differential equation: 

dX t = -aX t dt + adW t H + dL t A 

X t >b for all t > 0, I (1.1) 
X = x, J 

where a € (—00,00), a € (0, 00), x € [6,00) and W H = {W^)t>o is a one-dimensional standard 
fractional Brownian motion on A with Hurst index H G (0, 1). Here, the process L = (Lt)t>o 
is defined to be the minimal, non-decreasing and non-negative process with initial value Lq = 0, 
which increases only when X hits the boundary b, so that 



( I(X t > b)dL t = 0, (1.2) 

JfO.ool 



'[0,oo) 

where !(•) is the indicator function. It can be shown that (see, e.g., |12l 130] ) the process L has an 
explicit expression as 

L t = max < 0, sup I —x + a X u du — aW^ 1 ) > , (1-3) 
[ se[0,t] \ Jo J J 

where Lt represents the (cumulative) local time of X at the boundary b. 

It is often the case that the reflecting barrier b is assumed to be zero in applications to queueing 
system, storage model, engineering, finance, etc. This is mainly due to the physical restriction of 
the state processes such as queue-length, inventory level, content process, stock prices and interest 
rates, which take non-negative values. We refer the reader to [12J and [30j for more details on 
reflected stochastic processes and their wide applications. 

Our interest lies in the statistical inference for the RFOU process (jl.ip . More specifically, our 
aim is to estimate the unknown drift parameter a € (—00, 00) in (II. lh based on observation of 
the state process {Xt}t>o- Recently, the authors of [6 J studied the maximum likelihood estimator 
(MLE) for the model when a G (0, 00) and H = 1/2 (i.e., the standard Brownian motion case), 
and established several important properties. The MLE ax of a, based on the process {Xt} up to 
a previously determined fixed time T, is given by 

bL T -^X t dX t 

- == exM ■ (L4) 
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The MLE 2t satisfies strong consistency and asymptotic normality as T —> oo. However, this 
estimator is biased and its mean squared error (MSE) depends on the unknown parameter to be 
estimated. We note that exact estimates for the bias and the MSE of the estimator olt are not 
available. As a remedy for this, a sequential estimation plan {r{h),a T ^) was proposed in [15]. It 
is assumed in [15] that the parameter a £ (—00,00) and the process {X{\ is observed until the 
observed Fisher information of the process exceeds a predetermined level of precision h (see also 
[7]). More precisely, {Xt} is observed over the random time interval [0, r(h)] where the stopping 
time r(h) is defined as 

r(h) : = inf U > : J X 2 s ds > h\ , < h< 00, (1.5) 
and the J 7 ^^ -measurable function ct T {h) defined by 



«r(/i) 



r(h) 

bL T (h) — I X s dX s 







(1.6) 



is a sequential estimator. Then the sequential estimation plan (r(h),a T ^) has shown to satisfy 
the following properties (cf. Chapter 17.5 of [19J or Chapter 5.2 of [3]): (a) it is unbiased; (b) the 
plan is closed, i.e., the time of the observation r(h) is finite with probability 1; (c) its MSE is a 
constant that does not depend on the parameter to be estimated; (d) not only it provides consistent 
estimation plan but also a T ^ is exactly normally distributed, which makes it possible to construct 
an exact confidence interval for the parameter a. 

Our main results are concerned with extending the aforementioned estimators (jl.4p . (jl.6p and 
their respective statistical properties to the case when the system is driven by a fractional Brownian 
motion with H € (0, 1/2) U (1/2, 1). We establish uniform exponential moment estimates of the 
RFOU process, which, in conjunction with certain integral representations and the fundamental 
martingales of fractional Brownian motions, leads to two types of fractional Girsanov formulas. 
Then, we obtain the standard MLE and prove its strong consistency and asymptotic normality. 
Furthermore, we derive the explicit expression for the sequential MLE and show that it is unbiased, 
uniformly normally distributed (over the entire parameter space which is the real line) , and efficient 
in the mean square error sense. 

The rest of the paper is organized as follows. In Section we provide a brief discussion 
and derivation about how the RFOU model can naturally arise in the applications. Section [3] is 
devoted to preliminaries on fractional Brownian motion and fractional calculus that are necessary 
in our analysis. In Section [H we prove an important transformation result between a fractional 
Brownian motion and a standard Brownian motion, which is crucial in proving the equivalence 
of the two derived Girsanov formulas in the later section. We then investigate in Section [5] the 
standard maximum likelihood estimator of the drift parameter and prove its strong consistency and 
asymptotic normality. In Section we derive the explicit formulas for the sequential maximum 
likelihood estimator and its mean squared error by assuming the process is observed until a certain 
information reaches a specified precision level. The sequential estimator is shown to be unbiased, 
uniformly normally distributed, and efficient in the mean square error sense. 

2 Motivation of the RFOU model 

Here we provide some details about how the RFOU model could arise in the applications. Firstly, 
in the context of financial time series modelling, the RFOU processes can be used to describe the 
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spot foreign exchange rate processes, the domestic interest rate processes, and even some asset 
price processes in a regulated financial market system (cf . H] ) with long range dependence and 
heavy tails stylized facts, which seem to be common to a wide variety of markets, instruments 
and periods [U [31]. Secondly, in engineering applications to queueing and storage systems, the 
RFOU model can play as the key approximating process for systems with reneging or balking 
customers (cf. |27l \28\ [29] and the references therein), subject to their inter-arrival and/or service 
time processes exhibiting long range dependence characteristics in the traffic data. In such cases, 
the drift parameter a carries the physical meaning of customers' reneging (or, balking) rate from 
the system. We provide a more detailed justification below. 

Consider a single-server, single-class queueing model under heavy traffic subject to their reneging 
behaviors. More precisely, we shall consider a sequence of single-server queueing systems indexed by 
N = 1, 2, . . ., and assume that the sequence of interarrival times {T^ — T^]}^ =1 are stationary 

with E(T^ N) - Tffjj) = l/A (iV) such that A (JV) -> A G (0, oo) as N -> oo. Also, we assume an 
invariance principle holds: 

where W H is a standard fractional Brownian motion with Hurst parameter H and EW H (t) 2 = t 2H 
and some scaling constant a > 0. For the sake of simplicity, we will assume A = 1. Then, under 
a mild assumption on Var(T5 — Ti_}) < K < oo, one gets the following functional central limit 
theorem for the arrival process A^ N \t) = XmLi <t),t> 0: 

where W H is a standard fractional Brownian motion (cf. Theorem 2.1 of [14]). 

With the arrival process {A^ N \t) : t > 0}, consider the queueing system operating with a 
constant service rate fj,^ > 0, if there are customers in the buffer, and otherwise the server 
becomes idle. Moreover, customers faced with long waiting times will abandon the system before 
receiving service; let R( N ^ (t) be the total number of customers who abandoned the system during 
the time interval [0, t]. Incorporating these conditions will yield the following equation on the queue 
length process (Q {N \t) : t > 0): 

QW(t) = Q^\0) + A( N \t)-R^(t)- [\ W I [Q{N){s)>0] d S 

J 

= QW(0) + #)(t) - jf* a^Q^\s)ds - + f\^I [Q(N)(s)=0] ds, 

where a^' > is a proportionality factor related with the customers' abandonment rate from the 
system. We assume that the abandonment factor is of 0(N H ), in particular, a'^' /N H — > a 
as — > oo for some constant a > 0. 
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In view of the functional central limit theorem scaling used in (|2.ip . we have that 



Q w ( Nt ) _ QW(0) AW(Nt) l f m a( m Q( m (c)dc 



o 



Nt 



» {N) Nt 1 , {N) 
NH~ + W] q M 7 [QWW=o] ds 

" AT^ + N H J^H Jq {S)ClS 

NH + N w Jo M J [QW( S )=o] ds - 

Lastly, we impose a "heavy traffic" assumption implying that the system processing capacity is 
balanced with the system load, that is, the "drift" term (A^' — fi^ N ') converges to zero, at a 
certain rate closely related with the scaling in (|2.ip : N 1 " 1 ^ (A^ — — > as N — > oo. With 

the initial condition — > x G [0, oo) and owing to the Lipschitz continuity property of the 
generalized Skorokhod (reflection) map [29\ [23] together with the continuous- mapping theorem, 

we finally obtain the weak convergence of the scaled queue length process {^-jpr^- '■ t > 0}^ =0 to 
the RFOU process with the reflecting boundary given by zero. We note that similar derivations are 
possible for the offered waiting time process, in the context of state-dependent admission control 
setup with customers' impatient behaviors (see [16] and the references therein). 



3 Preliminaries on fractional calculus and fractional Brownian 
motion 

3.1 Fractional calculus 

In this section, we recall some basic results from fractional calculus. See [26J for more details. Let 
a,i £ 1 with a < b and let a > 0. (The symbol a in this section should not be confused with 
the parameter of the RFOU process.) The left-sided and right-sided fractional Riemann-Liouville 
integrals of / G ^([a, b\) of order a are defined for almost all t G (a, b) by 

1 r* 



i: + f(t) = — / (t- s r^f(s)ds, 

r (a) J a 

and 

ib-f (t) = ^rr- T ( s - *) a_1 f ^ ds > 
r(a) j t 

respectively, where (— l)~ a = e~ t7Ta and V (a) = / r a ~ 1 e~ r dr is the Euler gamma function. 

Jo 

Let I%+(LP) (resp. I^_{L P )) be the image of L p ([a,b]) by the operator J° + (resp. /„+). 

If / G I" + (L P ) (resp. / G I£_(L P ) ) and a G (0,1), then the left and right-sided fractional 
derivatives are defined as 



D« + f (t) = — L- ( _ZW + a f /W-ZWd, ) (3.i) 
r(i-a) V (*-<*) A (t-s) a+1 1 
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and 

n? f (+\ — , , , , 

r(l-a) V(6-*) a Jt (s-t) a+1 

for almost all t € (a, b). 

Let C a ([a,b]) denote the space of a-H61der continuous functions of order a on the interval 
[a, b]. When ap > 1, then we have J° + (L P ) £ C a p([a,b]). On the other hand, if /? > a, then 
C0([a, 6]) C for all p > 1. 

The following inversion formulas hold: 

D% + (I2 + f) = f, feL 1 

i° + (D° + f) = f, 

D^D^f) = D^ff, f e /^(L 1 ), a + /? < 1. 

Similar inversion formulas hold for the operators /£_ and .D°_ as well. 
We also have the following integration by parts formula. 

Proposition 3.1. If f £ I2 + {L p ),g e I£_{L q ) and± + \ = l, we have 

f {D a a+ f){s)g(s)ds = f f(s)(DZ_g)(s)ds. 

J a J a 

The following proposition indicates the relationship between Young's integral and Lebesgue 
integral. 

Proposition 3.2. Suppose that f £ C x (a,b) and g € C /J '(a,b) with A + \i > 1. Let A > a and 
/i > I — a. Then the Riemann-Stieltjes integral fdg exists and it can be expressed as 



f fdg = (-I)" f D% + f (t) Dlz a g b - (t) dt . 

J a J a 



(3.3) 



where g b - (t) = g (t) - g(b). 



3.2 Fractional Brownian motion 

The main reference for this section is [2j Chapter 2]. The fractional Brownian motion {Wjr ,t > 0} 
is a Gaussian process with zero mean and covariance function 

E(W t H W?) = \(t 2H + s 2H -\t- s\ 2H ) =: R H (t,s). 

Denote by £ the set of step functions on [0, T]. Let % be the Hilbert space defined as the closure 
of £ with respect to the scalar product (I[o,i]> I[o,s]) = R(t,s). Then the mapping L Q t -i — > Wf 1 can 
be extended to be an isometry ip — > W (ip) between H and the Gaussian space generated by W H . 
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The covariance kernel Rjj(t,s) has the following integral representation 

f(As 

R H (t,s)= / K H (t,r)K H (s,r)dr, 
Jo 



where 

r t „ i „ i 1 1 „ /"* 

K H (t,s) = b H 



{lf-\(t - s) H - 1 * -(H- \)s k 2~ H / " u H ~l (u - s) H -hu 

S A 



1 

, where T(x) = / s x ~ e~ s ds for x > is the Gamma 

T(H + i)T(2-2H) J Jo 

function. 

In particular, when H € (|, 1), Kjj(t, s) can be all written as 

K H (t,s) = C H sz~ H (u- s) H ~2 U H -^du, 

J s 

with C H = [H(2H - l)//3(2 -2H,H- ±)]t = 6 H (# - |), where (3(x,y) = / ^(l - tf^dt for 
x > 0, y > is the Beta function. 

Let Kjj denote the operator on L 2 ([0, T]), which is an isomorphism from L 2 ([0, T]) onto 

i£ + *(L»([0,3m 

(Xtf/)(t) = / K H (t,s)f(s)ds. 
J o 

We can rewrite the action of Kh as (see jlOj ) 

(W)W = ^r(i?-i)/ 1 +s ^C~ l (^- i? /)(i), if h>\, 

(K H f)(t) = b H T(H+\)I™s 1 2- H ll; H { S H - 1 2m), X H< X ~. 
If / is absolutely continuous, for H < ^, we can write 

(K H f){t) = b H T{H + 1 -)ll + s H -hD^{sl- H f)(t). 
Consider the operator K* H from E to L 2 ([0, T]) defined as 

(K* H <p)(s) = K H (T, s)<p(s) + j\<p{r) - <p(s))^(r, s)dr. 

Noting that (K* h I^q q)( s ) = Ku(t, s)I[ q (s), K* H can be extended to be an isometry between % and 
L 2 ([0,T]), 

(f,g)n = {K* H f,K* H g) L 2 {[ Q )T]) , f,g £ H. 

As a consequence, the operator K* H provides an isometry between the Hilbert spaces H and 
L 2 ([0,T]). Hence the process 

W t = W H {(K* H )- l I m ), t>0 (3.4) 
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is a standard Brownian motion, and W has the following integral representation 

W t H = f K H (t,s)dW s . 
Jo 

Moreover, we have the relationship between the Wiener integrals with respect to fractional 
Brownian motion and its related Brownian motion as 

W H &) := [ T <p(s)dW s H = [ T (K* H <p)(t)dW t , V en. (3.5) 
Jo Jo 

4 Transformation between fractional Brownian motion and Brow- 
nian motion 

We shall establish a useful transformation result between fractional Brownian motion and Brownian 
motion. Let < s < t < T, and let 

k H {t, S ) = K- 1 S 1 2- H {t-s) 1 2- H , 

where k h = 2flT(| - H)T(H + \). 
Define 

Mf = f k H (t,s)dW s H . 
Jo 

Then using Equation (|3.5p . we may get that M is a Gaussian martingale with quadratic variation 

(M H ) t = \ H H 2 ~ 2H , where X H = ( ~7_ ' 2 - The martingale M H was introduced as the 

fundamental martingale for the fractional Brownian motion W H in order to get an Girsanov type 
of theorem for W H in [20] . 

By Levy Characterization Theorem, it is easy to verify that the process 

Bt = H(2H-i) r* a h _ Um ^ t ^ Q (41) 

Jo 

is a standard Brownian motion. 

On the other hand, we know that W H can be represented as an Ito integral with respect to the 
Brownian motion Wt = W H ((K^) -1 1[o,t]) ■ A natural question arises: does B coincide with W? 
The answer is positive and, in fact, we have the following theorem. Hereafter, we shall use C to 
denote a generic positive constant, which may vary from line to line. 

Theorem 4.1. Let the process B = {B{\t>o be defined as in (|4.ip . Then B is the standard 
Brownian motion related with W H , i.e., 

W t H = f K H (t,s)dB s . 
Jo 

Proof. Denote 

Z t = kJj 1 / s^~ H dW 11 , t > 0, 
Jo 
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and hence we have 



and 



W t H = K H I s n ~idZ. 



M 



rt 1 
f = / {t-r)" H dZ r , s > 0. 
J o 



An equivalent equation for the generalized Abel equation (14. 2\\ is given by 

Z t = — r i-. !\t - r) H ~ldM?. 

In fact, when H > 1/2, to solve (|4.2p for we observe that 



(t - s) H -2M^ds 



I (t-s) H -^ [ (s -r)^~ H dZ r ds 
Jo Jo 
ft ft 

/ / (t- s) H -i(s-r)2~ H dsdZ r 

Jo Jr 



P{H- l -*--H)Z t . 



Equation (|4.3p follows by the fact 



(t-s) H -2M?ds 



2H - 1 



(t-s) H "2dMf . 



When .ff < 1/2, we may get (j4.2|) from (14. 3p using the same strategy as for H > 1/2. 
Now we show that 



f s H -idZ s = t H ^Z t -(H- 1 -) f s H -lz s ds. 
Jo 1 Jo 



(4.2) 



(4.3) 



(4.4) 



When H > 1/2, (|4.4p is a direct consequence of integration by parts formula. When H < 1/2, 
notice that ^ 

/ s H ~^dZ s = t H ~^ Z t -lim ( s H -^ Z s - (H --) I r H -lz r dr). 
Jo S 4-Q V 2 is / 

Hence, it suffices to show that \\ms H ~^ Z s = 0, which implies that / s H ^^Z s ds is well-defined 

slO J Q 

and therefore (|4.4p holds. In fact, we have 



C 



A- H dw* 



C (\s^- H W s H \ + J \r-5- H W?\dA <Cs^ H \\W Hl 



where the norm II • lira , is defined in (15.21). and \\W 



oo,s 

Yuii\s h -^ZJ < CHW^Hooilim.*" - 0. 



oo,i- Consequently, we obtain 

H 



slO 



40 
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From integration by parts, Fubini's theorem and the fact dM^ 



C H 



r2 H dB r , we have 



h(2h - r 



W? = K H 



Jo 



dZ, 



K H 



1 



t H —*Z t -(H--) / s H ~^Z s ds 
2 Jo 



1 



r(ff + i)r(§-ir 



■ [\t-r] 

Jo 



H ~^dM? 



T{H + \)T{l-H) 



s H ~2 I (s -r) H -2dM?ds 



H(2H - 1) 



r(ff + i)r(§-fr 



; ?r \ r / f s H ~Us -r) H -^dsA~ H dB r 

= b H f K H (t,r)dB r . 
Jo 



□ 



5 Parameter estimation for the RFOU 

We begin by obtaining a uniform moment estimate for the RFOU process. Then we derive two 
kinds of Girsanov formula and prove their equivalence when H > 1/2, which leads to the maximum 
likelihood estimator of the drift parameter. 

5.1 Estimate on the RFOU process 

For any (3 € (0, 1), we denote by C@(0, T) the space of /3-H61der continuous functions on the interval 
[0, T\. We will make use of the notations 

II II \%r ~ X d\ /r 1 \ 

\\x\\ pT = sup — — —s-, (5.1) 

0<8<r<T \r - V\ p 



and 



Mloo,r= sup |av|. (5.2) 

0<r<T 



Lemma 5.1. Let X t be the solution to the following Skorohod equation 

dX t = -aX t dt + adW t H + dL t 

subject to X(0) = x > 0, where L is the minimal non- decreasing process. Then X is (H-) Holder 
continuous. There exists Ao > depending on T, H and e such that 

Eexp(Ao||-X'||^ 0)T ) < oo and Eexp(Ao||-X"||^-_ eT ) < oo. (5.3) 



10 



Proof. Note that from (jl.3p . we have 



X t = x - a [ X s ds + aW t H + sup {-x + a [ X u du - aW^} V 0. 

Jo 0<s<t Jo 



Assume a > 0. We have that 

X t < x - a [ X s ds + crW t H + a [ X s ds + sup {-x - aW^ 1 } V 

Jo JO 0<s<t 

< 2x + 2a\\W H \\ 00tT . 
For general a £ M, we have 

\X t \<2x + 2\a\ f \X s \ds + 2a\\W H \\ 00>T . 
Jo 

By Gronwall's inequality, we have 

\\X\\ 00)T <2(x + a)e 2 ^ T \\W H \\ 00)T 

Then, the first inequality in (|5.3[) follows from the Fernique's Theorem, which claims the exponential 
integrability of the square of a seminorm for a Gaussian process (see 

For the second inequality, noting that the function f(x) = xV0:= max{x, 0} is Lipschitz, and 
that for a general continuous function g, 

sup {g s } V - sup {g s } V < ( sup {g s } - sup {g s } ) V 

0<s<t+r 0<s<t \t<s<t+r 0<s<t J 

< sup {g s } - 9t 

t<s<t+r 

< sup {g s - g u }, 

t<s,u<t+r 

we have 

\X r - X s \ < 2|a|||^:|| 00 (r - s) + 2a\\W \\ H -e, T (r - s) H ~ e . 



So we get that 



\X\\ h _ 6 t < 2|a|||X|| 00 , T T 1+e - H + 2a\\W\\H-e,T. 



Applying Fernique's Theorem once more, we have lEe XllWll »-^T < oo for some A > 0. Combining 
with the first inequality for HAH^t, we have the second inequality. □ 

5.2 MLE based on the Girsanov transform I 

Here is the Girsanov theorem for the fractional Brownian motion, the version from |22j . 

Lemma 5.2. Let W t H = W t H + / u s ds , where {u t ,t G [0,T]} is an adapted process with respect 

Jo 

to the filtration J :W generated by W H with integrable trajectories. Denote 

6r = exp ^- £ (k^ jf u r dr^j (s)dW s (k^ u r dr^j (s)dsj , 

where W is the Wiener Process defined by (TO^. Assume that 



11 



(i) I u s ds £ I Q+ 2 (L ([0,T])) ; almost surely. 



(ii) E(fr) = 1. 

Fhen the shifted 
the new probability P defined by dP/dP = £j\ 



Then the shifted process W H is an -fractional Brownian motion with Hurst parameter H under 



The following is the Girsanov type theorem for RFOU. 
Theorem 5.3. Let 



£t = exp ( J I if^ 1 J ^X r drj (s)dW s -^J iK^ 1 J -X r dr ) (s)ds 



T 



a 



Then {X s : < s < T} is a reflected fractional Brownian motion under the new probability P 
defined by dP/dP = 



Proof. If we can show that E{fyc) = 1, then by Lemma E3 W t H = Wf - - / X s ds is a fractional 

_ °~ Jo 

Brownian motion under P. Hence X% = Xq + aWt + Lt is a reflected fractional Brownian motion 
under P. We now show that for At small enough, we have 



rt+At / /•• \ 2 \ 

E(>xp( / IK' 1 J X X r dr) (s)dsj<oo. 



(5.4) 



We shall prove this in two cases H < \ and H > \, respectively. 



Case H < 2- Observe that 



K 



H 



X r dr (s) 



s H -ht H s^ H X s 



l 0+ 



Cs H - L * 



(s - r)-^- H r^' H X r dr 



< C(l + 1x1 + 11X1100^). 



Then when At small enough, (|5.4p follows by Lemma 1571 
Case H > 2- We have that 



K H l J Xrdr) (s) 
s^ H X s + (H - \)s H -h ( Xs [ 8 *~*_~^~* dr + 



o (s — r) 1 - 



X s — X r 
o (s — r)^ +H 



r 2 H dr 



Noting that 



S \-H _ r \-H 



1+H 



dr = Cs 1 - 2 ", 



o (s — r) 
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X s -X r \ < \\X\\ H -e,T( s -r) 



H- 



< r < s < T, 



and 

we have, if we choose e < ^, 

K^ 1 ^ X r dr) (s) < C ^ H s H ^\\X\\ H ^ T + s H -^\X\\ H . e , T J\s -r) H ^(s -r)- l 2- H r^ 

< C\\X\\ H _ e , T {s^ + s H ^) 

< C\\X\\n-e,T- 

Then when At small enough, (|5.4p follows from Lemma 15.11 

T iT 

Now choose n big enough such that — < At, and let t, = — , i = 1, ••• ,n. Denoting 

n n 

(Kh 1 Jo % x rdr) (s) by f(s), we obtain 



exp f(s)dW s -±£ f 2 (s)ds^j 



ie I it: 

IE 



(/ ftn-1 1 ftn-1 

exp (J f( s )dW a --j f{s)ds 



We have that iB 



xp ( f T f(s)dW s -\\ T f\s)ds) 
exp ( / T /( S )dW s - i / T / 2 ( S )cfe 

\ *n— 1 *n-l 



< 1 a.s., since 



1 



exp I / f(s)dW s / f 2 (s)ds is a positive local martingale and hence a supermartingale. 



On the other hand, since iB 



1 

exp - 



/ 2 (s)ds 



tn-l 



< oo, then by Novikov Criterion, we know 



that iE 



exp 



/(«)dw s - i r f 2 (s)d S 

tn — l •* tn—1 J 



1. This implies that 



£7 



exp 



f(s)dW s 



T 



f\s)ds 



tn-l 



1 a.s. 



We use the above procedure n times, and we may get that 1E{^t) = 1- 



□ 



Denote by (C([0, T]),Bt) a measurable space of continuous functions {x s ,0 < s < T} with 
initial value xo = 0. Let X R = {X R (t)}t>o be the reflected fractional Brownian motion, that is, 
the model (jl.ip with a = 0. Let P a and P R denote the probability measures on (C([0, T]), Br) 
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induced by X and X R , respectively. Similar to the proof of [18, Theorem 7.1], we can show that 
P a ~ P R , and 

dP a . 1 / / x /„ i /'a.. . \ 1 f r f „ 1 /"a 



= exp f - jf jf ^X r dr^ (s)(iW s + i jf (k^ j ^X r dr^j (s)ds 



Let 4>(t) = / X s ds, and we have 

r (k^v) (s)dw s = f {{k^k^) ( S )dw s H . 

Jo Jo 
The integral above is in the sense of Skorohod. Formally, we may represent a Skorohod integral as 
the sum of a Stratonovich integral and a correction term involving Malliavin derivative. 

Noting that dWf = ±dX s + ^X s ds - ±dL s , we have 

rp 

[ {Kg 1 <£) (s)dW s 

Jo 

= [ {(K* H )- 1 K H 1 cl ) )(s)odW s H + E T 
Jo 

= -f ((K* H )~ 1 K] i 1 cl ) )(s)(dX s + aX s ds-dL s ). 
a Jo 

where 

E T = F \(K* H )- 1 K H 1 <p)(s)dW s H - [ T ((K* H y 1 K H 1 ^(s)odW s H , 
Jo Jo 

the symbol o meaning Stratonovich product, and if (K^^K' 1 ^ G B 1 ' 2 (|^|) then 

r T r T 

E T = H(2H - 1) / D a ((K^-^-K- 1 ^) (t)dsdt (see, for instance, [STJ Proposition 5.2.3]) 
JO JO 

Note that when H > ^, 

F {{K^K^cf) (s)X s ds = F {(K^K^cf) ( s )</>(ds) = F {Kj^cf) (s) 2 ds, 
Jo Jo Jo 

where the last step follows [2j Lemma 2.1.9]. We then obtain that 

^ = exp (-£ [ (m)-^ f o X r dr) { s)dX s ~\Q 2 [ jf X r dr) % (s)ds 

+^ £ ((K^Kh 1 J X r dr^j (s)dL s + E T ^j . 

dP a 

The MLE for a given the realization X is the value which maximizes ( see ! e -g-> |23j for the 

definition of MLE). So the MLE a T is 

/ r ((K^K- 1 J Q X r dr) (s)dL s - J T ((K^K^ £ X r dr) (s)dX s + E T 

OLT '■= Tfi 9 • 

ff{K?&Xrdr) 2 (s)d8 

However, St is not a practical statistical estimator, since it is infeasible to evaluate the quantity 
Et- We shall obtain a more practical MLE using the Girsanov theorem obtained via a fundamental 
martingale approach in the next subsection. 
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5.3 MLE based on the Girsanov transform II 



We refer the reader to [20] for the Girsanov transform for a fractional Brownian motion based on 
the fundamental martingale. Let < s < t < T, and recall (from Section H]) that 

k H (t, s) = K -£ a *- H {t - s^' H , k h = 2HT(- - H)T(H + -); 

M t H = [ k H (t,s)dW*. 
Jo 

The quadratic variation of a martingale M H is (M H ) t = Xjj t 2 ~ 2H , where A# = rd-m • 

Let 

X t = f k H (t,s)dX s , L t = f k H (t,s)dL s , and Xt = ttttwv / k H (t,s)X s ds. (5.5) 
Jo Jo d \ M it Jo 

Then we have 

dX t = -axtd(M H ) t + adM t H + dL t . (5.6) 
Lemma 5.4. There exists a constant C > depending on H and e such that 

f X 2 s d(M H ) s <C(T 2 - 2H \\X\\ 2 ^ T + T 2 -z\\X\\ 2 H _ £tT ) . 
Jo 

Proof. We begin by noting that 



Let Y t = S 2- H X s ds, then we have ||y||oo,T < CWX^^T^ 11 . 
Jo 



1 

Case 1: H > -. We have that 



Jo 



S )2~ H S 2- H X,ds 



\t-s) 1 2- H d(Y s -Y t ) 



o 



lim(* - S )i" H (y s - Y t ) + tl~ H Y t + (H-h [ (Y t - Y s )(t 
sft 2 J 

t\- H Y t + {H-\) j\Yt-Y s [ 

t l 2- H Y t + (H-h f{Y t - Y t ^ s )s-^- H c 
1 Jo 



,)(t - s)~^~ H ds 



Also, observe that 



^ [\Y t -Y t „ s )s^- H ds 
at Jo 

lim ± ^ £ {Y t+£ - Y t+ ^ s )s-^- H ds ~ J Q ( Y t- Y t _ s )s-h- H ds^ 



--: A(t)+B(t), 
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where 



A(t) = lim X - ^ + \Y t+£ - Y t+£ „ s )s^ H ds - j\v t+e - Y t+£ ^ S ) S -^ H ds^j , 



and 



B(t) = lim X - (J (Y t+£ - Y t+£ ^ s )s-^- H ds - J (Y t - Y t _ s ) S -h- H ds^ . 

It is clear that A(t) = t~i~ H Y t , when t > since Y t is continuous in t. For B{t), we have 

f(Yl-YU)s-l- H d S 
Jo 

f (th- H x t -(t-s)^ H x t _ s ) s-^ds 

f (X t - X t ^ s ) s-^ H ds + jf X t _ s (ti~ B -(t- s)h~ H ) s~i- H d ; 



B(t) 



Observe that 



0< f (X t -X t - s )s-^- H ds< [ \\X\\ H ^ T s H - 
Jo Jo 



h-H^_nx\\ H - 



■S 2 M ds 



■e,T ,i- 



and also that 



f t ^ H ~(t- s^- H s -^- H ds = t 1 - 211 r i - a - S ) 

Jo Jo 



l-2e 



_l_rr , 
S 2 " < OO 



since lim 

■40 



th~ H - (t-s)^- H 



— - — H 

s 2 < oo. Hence 



|B(t)| < C (HXllH-e.Tt 1 ^^ + ||Jr||oo,T* 1 - 2H ) • 



Therefore, 



d f\ 



s)\~ H S \- H X s ds 



C 



t-^- H Y t + ti- H ti- H X t + B{t) 



hA-h 



and for Xt we have 
Hence, we conclude that 



< c (\\x\UTt 1 - 211 + \\x\\ H _ etT t l - H -*) 
Ixtl^ciwxw^T + wxWn^Tt 11 -*). 



f X ld(M») s = \ H \2-2H) f xy- 2 ^<C(||X||^ 2 - 2ff + ||X||^ £jT t 2 -) 
Jo Jo 



Case 2: H < -. Notice that 
2 



d_ < > 

dt 



= r^-ioi^-V-^Xf), 
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and also 



I> + (s*- a X a )(t) 



r(i-F) 



(t - S yi- H S2~ H X s ds 



< C\\X\\^ T t 



1-2H 



Therefore, similar to the results in Case 1, we have 

\Xt\ < C||X||oo,T. 

Hence, we get 

l\ 2 s d{M H ) s = \ H \2 - 2H) l\y- 2H ds < C\\X\\l tT i 
Jo Jo 

Remark 5.5. From the proof the above lemma, we have that 



2-2H 



□ 



• ifH> -, 



d_ '■' 

dt 



J\t - s)^- H sh- H X s ds = t l i- H Y t + (H-±)(A(t) + B(t)) 



t l - ZH X t + {H --)B{t) 



t^ H X t + {H-\) jf* [t*- H X t - (t - sf>- H X t _ s ) s-^ H ds 



H-i, I. 



r(--H)D 0+ *(u-2- H x u )(t) 



j t J\t - s)^ H sl- H x sds = r(| - H)il; H (J-»x u )(t). 

If we define Dq + as Iq® for a < 0, then for all H € (0, 1), xt can be represented as 

ni- h)\ h 2H _ h-i , 



Xt 



2-2H 



K^t 2 "-'D» + *(u2-"X u )(t). 



Remark 5.6. When H > 



-t H -±D^{^- H v'){t), 



and when H < -, 



.1 



H.J\ 



b H T(H + 



17 



By a direct computation, we can get 

rT rT 



[ X tdM t H = [ (K H l ! X u du)(s)dB s , 
Jo Jo Jo 



where B is the Wiener process defined in fl^. By Theorem \4-l\ we know that B coincides with 
W, the Wiener process related with the fractional Brownian motion W H . This implies that, the 
two types of Girsanov theorems, Theorem 15,31 and Theorem 5. 1 below, coincide with each other. 

Theorem 5.7. Let 

rT „,2 rT 



Trr = w(^j XsdMf-^J x 2 s d{M H ) 



Then {X s : < s < T} is a reflected fractional Brownian motion under the new probability P 
defined by dP/dP = rjr- 



Proof. It suffices to prove that JE(t]t) = 1, which is a consequence of Lemma 15.41 and Lemma 
15.11 for sufficiently small T. For general T, we just split the interval [0, T] into sufficiently small 
subintervals, and employ the strategy used at the end the proof of Theorem 15.71 □ 



Let P a and P R denote the probability measures on (C([0, T]),Br) induced by X and X R 
respectively, where X R is the reflected fractional Brownian motions (i.e. a = 0). Similar to the 
proof of [18J Theorem 7.1], we can show that P a ~ P R , and 



dP a 
dP^ 



(rT 2 rT 

-^J^ XsdMf + ^j^ X 2 s d{M H ) 1 



exp 



a 



XsdX s 



a 

2^ 



xid(M% + 



a 



XsdL s . 



Hence the MLE for a is given by 



ax '■- 



Jo XsdX s + f T x s dL s 



x 2 s d(M«) s 

Before we prove the consistency of ay, we need the following lemma. 

Lemma 5.8. Let {Y t ,t > 0} be the fractional Ornstein-Uhlenbeck (FOU) process satisfying 

dY t = -aY t dt + adW t H 

with initial value Y$ = x > 0. Then we have Xt > \Yt\ a.s. for t > 0. Furthermore, we also have 



lim i ( \Y t \dt = lE a [ e^-^cWf + 
^oo T J [ J_ 00 

1 f T 

hm — / \Y t \dt = oo, 



x 



when a > 0, 
when a < 0. 



Proof. Since d(X t - Y t ) = -a(X t - Y t )dt + dL t with X -Y = 0, we have X t > Y t a.s. for all t. 
Let Zt be the FOU process satisfying dZt = —aZtdt — adW^ with initial value Zq = —x. Then 
we have Z t = —Y t . On the other hand, d(X t + Z t ) = — a(X t + Z t )dt + dL t with Xq + Zq = 0, and 
hence — X t < Z t = —Y t . So we have \Y t \ < X t a.s. for all t. 
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Note that Yt has the following expression, 

pt r-t fO 

Y t = x + a / er a{t - s Uw^ = x + a / e~ a{t - s) dW s H - e~ at a / e as dW s H . 

JO J— oo J— oo 

When q > 0, the process |lj = a J e~ a ^~ s ^ dW^ ,t > o| is Gaussian, stationary and ergodic 
(see [5]), and by ergodic theorem, 



which implies that 



lim - / |i^ + a;|dt = ^(|yo + a;|), 

T^oo i J 



lim 1 / |y t |dt = ^(|F + a;|) > 0. 



When a = 0,Y t = aW t H ', and let y t = a [ e s dWf . Note that e~V = erW/* - / e~ s y s ds, and 

Jo Jo 



Jim c-'ijkII IL'[I ae s dWf + .i j r. 



(5.7) 



Hence we have 



> 



1 /" , , 
liminf— / Yf di 

T^oo T7 ' " 
1 ' T " 



liminf — / / e s v s dsdt 



JO 



1 / ,T 

lim - / e l y t dt 



OO, 



1 f T f f - 

where the last equality holds because of (|5.7p and the fact that liminf — / e s y s dsdt = oo 



since / e s y s ds behaves like xt when i goes to nifty. 
Jo 

When a < 0, let j3 = —a > 0, and then 



ft fOO poo 

Y t = x + a e^dWf = x + a e^dW, 11 - a / eP^dW?. 
Jo Jo Jt 



Similarly as in [8], we can show that | Yt = o J eP^ s ^dW^ , t > Oj is stationary and ergodic. 
Note that now Yt = x + e^*Yo — Y and hence 



f T \Y t \dt > \Y \ lim i f T e^dt- lim i C \x-Y t 



liminf- / l^ldi > 

> oo a.s 



\dt 



□ 



Theorem 5.9. For H € (0, 1), the MLE cut is strongly consistent. 
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Proof. To prove lim ay = a a.s., it suffices to show that 

T— too 



Since / XsdM^ is a martingale, by Lepingle's law of large numbers (see [18, Lemma 17.4]), it 
Jo 



suffices to show that 



From Remark 15.51 we know that 



lim 

T->oo 



X 2 s (s)d(M H ) s = oo. 



(5.8) 



H-- --H 

Recall that D 0+ 2 actually means I 2 + when H < ^. 



Let c = 1 - 2H and d = AH - 2. Then c + d = 2H - 1 and c + § = 0. Notice that 



JO 



C / s 
'o 

c 



2H-1 



o 



D 0+ 5 n2-^X u ) (s)ds 

2 



s | si ( Do + l2 u2- H X u ) ( s ) ) ds 



> C 
= C 

= c 



1+c , 



r 



S - ' 2 



D 0+ 2 u2- H X u {s)ds 



[ (D^Jl](s)s 1 2- H X,ds 



T~ Jo 
1 

rpl-H 



(T - s)^- H s^~ H X s ds 







where the only inequality above follows by Jensen's inequality or Holder's inequality. 
1 



If H > — , we have 



-H 



T 



(T - s)^- H sz- H X s ds 



1 2 


h- > 




[T H Jo 



> c 



T H Jo 



X*ds 



Therefore, to prove (|5.8p . it suffices to show that 



1 f T 
lim — / X s ds > 0, 
T J 



(5.9) 



and it is an immediate consequence of Lemma 15.8 

If H < — , we can show that 
~ 2' 



1 f i i 

liminf^ / (T - s)^ H s^ H X s ds > 0, 



(5.10) 
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which is a sufficient condition for inequality (|5.8p . 
In fact, we have 

fT 



(T - s)^- H s^~ H X s ds 
rT / r T , , 

S 2~ H X,ds 







= C [ (( (T-r)-^ H dr 

JO \Js J 

= C £ ^ jf" s^- H X s ds^j (T - r)^- H dr. 

For the term / s2~ H X«ds, we have 
Jo 

i r i i /" r i i /" r i /" r 

liminf — / s~ H X s ds > liminf — / sz~ H X s ds > liminf — / X s ds = lim — / X s ds 

r— >oo r Jq r— »oo r > — >oo r J i r— >oo r Jq 



where a is a positive number or a = oo by Lemma 15.81 

When a is finite, for any e > 0, there exists N > 0, such that — / s^~ H X s ds > a — e for all 

r 



r > N. 

Therefore 



o 



liminf— = / (T — s)^~ H s^~ H X s ds 

T^oo T2~ H Jo 

= liminf— ^ — / ( f s^~ H X s ds) (T - r)"^~ H dr 

T l- H JN \Jo J 

> liminf — g / (a — e)r(T — r)~^~ H dr 

T^oo Tl~ H Jn 

=P(l,±-H)(a-e), 

which implies (|5.10p . 

When a = oo, similarly we can show ()5.10p holds with oo as the value of limit on the left hand 
side. □ 

Theorem 5.10. For H £ (0,1), the MLE a? of a admits the asymptotic normality, i.e., 



aT ^J j\^d(M H ) s A N(0,1) asT^oo. 



Proof. This is an immediate consequence of the Central Limit Theorem for martingales (see, for 
instance, [231 Theorem B.10].) □ 



6 Sequential MLE 

Recall the strong consistency and asymptotic normality of the MLE ax, established in Theorems 
15. 9115.101 Such results are valuable from statistical analysis viewpoint of applications, however, 
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generally speaking, the MLE is a biased estimator and its mean squared error (MSE) depends 
on the parameter to be estimated (Theorem 15. 10H . Moreover, the classical Cramer- Rao lower 
bound may not be attained for this MLE. To overcome such limitations, we consider the sequential 
estimation plan and verify that the proposed plan is significantly helpful both in asymptotic and 
non-asymptotic short time observation. 

In contrast to the MLE, the proposed sequential maximum likelihood estimator (SMLE) is 
unbiased, exactly normally distributed (on the finite time observation), and its MSE has an explicit, 
simple expression that does not depend on the parameter to be estimated (see Theorem 16.11 below) . 
The SMLE is uniformly normally distributed over the entire parameter space which is the real line. 
Such results would be of ample use in applications to several areas such as engineering, financial 
and biological modeling where unknown parameter estimation is based on relatively shorter time 
observation. Furthermore, an analog of the Cramer-Rao lower bound is proved and the SMLE is 
shown to be efficient among all unbiased estimation plans in the mean squared error sense (see 
Theorem 16.21 below) . 

Define the stopping time r H (h) as 

r H (h) := inf jt > : xld(M H ) s > /ij , < h< oo. (6.1) 
Then, J 7 ^,^ -measurable function a T HQ^ 



a T H (h) 



T H (h) fr H (h) 

XsdL s - / XsdX s 
o Jo 



(6.2) 



is a sequential estimator. From (|5.5p ^ (|5.6p . it can be seen that 



<*T«(h) =a ~\\ XsdMf. (6.3) 



o 



A proof of the next theorem follows along the similar lines of [19J by accommodating our basic 
model assumptions involving the reflection (i.e., state space) constraint and fractional Brownian 
noise with H > 1/2. In what follows, the index a in P and E emphasizes the fact that the 
distribution of the state process is being considered for the prescribed value a. 

Theorem 6.1. The sequential estimation plan (t h (h),a T H (m) , < h < oo, has the following 
properties: 



(i) P a (T H (h) < oo) = 1 for all a G (-oo, oo), 

(ii) K a (a T H f h \) = a for all a E (— oo, oo), 

(Hi) E a (a T H {h) - a) 2 = —, 

c o- 2 

(iv) a rH(h) - a = N(0, —), 

(v) a T HrfA — > a a.s. as h — > oo. 
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Proof. We first show that the J~*H(h\ -measurable random variable a T n^ is indeed the SMLE as 
follows. Consider an arbitrary stopping time r with respect to the filtration {J 7 i^}t>o generated by 
the process X. Similar to |18|, Theorem 7.1], the probability measures ^ x anc ^ x induced by 
the processes stopped at time r 

dXf = -OXfdt + adW t H + dL 9 t , X? > b for all t > 0, 1 
dXf = -aX?dt + adW t H + dLf, X? >b for alU > 0, J 

respectively, are equivalent and their Radon-Nikodym derivative is given by 



dKx 



\X 



T.fl 



exp |-i J\a - 6) X ?dM t H + ± J\<* 



(6.4) 



where is the natural filtration generated by {Xf : < t < r}, and 



d(M H ) 



t Jo 



k H (t,s)Xfds. 



Let X® be the reflected fractional Brownian motion (RFBM) satisfying dX® = adWp + dL t , 
t > and P°- be the measure induced by the RFBM X°. Then the log likelihood function £ T (a) is 
given by 



£ T (a):=a 2 log—^ = -a x fdX t -^- ( X f) 2 d(M H ) t + a / X t^ t . 

dr r,X JO 1 JO JO 

Then, by solving the equation 



d ( , d¥ a 



we obtain the SMLE given by 



ot T = Si -= o = a = — - — n • 0.5) 

f Q T (xf) 2 d(M H ) s fo (x?) 2 d(M») s 

Now, setting r = r (/t) in (|6.5p . we obtain for S r = a T H/ h \ the representation given by (|6.3|) . 



We shall use x instead of x a - To verify part (i), it suffices to show that lim L x 2 d(M H ) s = 

T— »oo 

a.s., which has been proved in Theorem 15.91 The claims in (ii),(iii) and (iv) follow from the fact 

that (Jq ^ XsdM^ 1 ,h > 0) is a standard Brownian motion indexed by h (see, for instance, \13\ 
Theorem 4.6]). Also, the result (v) is an immediate consequence of (iv). □ 

Next, we show that the proposed sequential estimation plan (t (h) , a T H (m) is efficient among 
all unbiased estimation plans in the mean squared error sense. To do this, we will assume that 
the reflecting boundary 6 = and prove an analog of the Cramer-Rao lower bound for arbitrary 
unbiased estimation plans. 
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Theorem 6.2. Assume 6 = 0. Let the sequential plan (r,a T (X)) be an arbitrary unbiased estima- 
tion plan for the RFOU process {Xt} with the parameter a € (—00, 00), namely, 



E a (a T (X)) = a for all a £ (— 00,00). 



Suppose also that < E^/J" x 2 d(M H ) s ] < 00. Then, 

Var Q (S T ) = Eo,[S r — a] 2 > 



a 



(6.6) 



(6.7) 



Remark 6.3. A sequential estimation plan (r,a T ) is said to be efficient in the MSE sense if for 
which (|6.7p becomes an equality for all a E (—00,00). Notice that since E a [a r H^ — a] 2 = 

as established in Theorem \6.1\ (Hi) and E a [/J" ^ Xsd{M H ) s ] = h by the definition ofr H (h), the 
sequential estimation plan (t h (h),a T H t^) is efficient in the MSE sense. 

Proof. Without loss of generality, we assume that a = 1. In view of the Radon-Nikodym derivative 
expression in (|6.4p with 6 = 0, differentiating both sides of (|6.6p with respect to a yields that 



E„ 





. L 







(cf. the proof of Theorem 7.22 in [19J). Then, since 



Xs dX s -a x 2 s d{M H ) s + Xs dL 



E 
0, 



1 



(6i 



Xs dX s + a I xtd(M H ) s - I Xs dL 



XsdM 1 



it follows that 



Xs dX s -a Xs d(M n ) s + Xs dL t 
Jo Jo 



(a T — a) 

Applying Cauchy-Schwarz inequality in (|6.9p . we obtain 
1 < E a [a T -a] 2 E a 

= E a [a T -a] 2 E a 
= E a [a T -a] 2 E a 



1. 



(6.9) 



XsdX s + a / xtd(M u ) t 
Jo 



H 



XsdL, 



XsdM s 
X 2 s d(M H ) 



(6.10) 



where the first equality follows from the state equation (jl.ip . Since < E Q [/ T x1d(M H ) s ] < 00, 
we can divide both sides of (|6.10p by this factor, and then the desired result follows. □ 
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